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Abstract — The  moments  produced  by  the  aerodynamic  con¬ 
trol  effectors  on  aircraft  are  typically  non-linear  functions  of 
displacement  with  a  region  of  linearity  about  zero  displacement. 
For  nominal  operation  with  no  locked  or  boating  effectors  and 
for  gentle  maneuvers,  affine  approximations  to  the  moment 
displacement  curves  are  sufficient  as  the  effectors  operate  pri¬ 
marily  in  the  linear  region.  However,  if  an  actuator  has  become 
locked  and/or  aggressive  maneuvering  is  required,  the  control 
effectors  may  require  large  debections  into  the  non-linear  regions 
of  the  moment-displacement  curve  in  order  to  maximize  the 
available  control  power.  As  we  move  away  from  the  linear 
region,  the  affine  approximation  becomes  inaccurate,  especially 
in  cases  where  the  slope  of  the  moment  curve  changes  sign.  In 
this  paper  we  will  brieby  review  a  novel  approach  that  was 
developed  at  AFRL  for  solving  the  control  allocation  problem 
where  the  control  moments  are  piecewise  linear  functions  of 
the  control  displacement.  The  advantages  of  this  approach  are 
that  it  avoids  cumbersome,  high-order  polynomial  curve  bts  of 
the  aerodynamic  data  and  thus  avoids  solution  of  a  non-linear 
programming  problem.  Stability  and  control  derivatives  are 
almost  always  stored  in  multi-dimensional  look-up  tables  where 
it  is  assumed  that  the  data  is  piecewise  linear.  The  approach 
that  is  presented  utilizes  this  piecewise  linear  assumption  for  the 
control  effector  moment  data  and  accounts  for  non-linearities  in 
the  moment-effector  relationships  and  position  constraints  on  the 
effectors  as  well.  The  control  allocation  problem  using  piecewise 
linear  functions  is  posed  as  a  mixed-integer  linear  program 
(MILP)  and  solved  using  a  freely  available  GNU-licensed  mixed- 
integer  linear  programming  code. 

We  will  demonstrate  two  applications  of  the  piecewise  linear 
control  allocation  (PLCA)  approach.  The  brst  application  to 
be  considered  is  to  use  the  PLCA  approach  in  the  inner- 
loop  control  law  of  a  re-usable  launch  vehicle  on  approach 
and  landing.  Body  axis  angular  rates  are  controlled  using  a 
dynamic  inversion  controller.  The  vehicle  will  be  subjected  to 
two  stuck  control  effectors,  and  recovery  of  the  vehicle  using 
only  control  effector  reconbguration  provided  by  the  control 
allocator  and  without  trajectory  reshaping  will  be  demonstrated. 
However,  the  downside  of  the  MILP  formulation  is  that  the 
optimization  problem  cannot  be  solved  fast  enough  for  real¬ 
time  implementation  on  the  current  generation  of  bight  control 
computers. 

The  second  application  that  will  be  demonstrated  is  constraint 
estimation  for  trajectory  reshaping  and  re- targeting.  In  order 
to  successfully  re- target  a  trajectory  when  an  aircraft  has 
experienced  degraded  performance  due  to  a  failure  or  damage 
to  the  vehicle,  the  effects  of  the  failure  or  damage  on  the  lift, 
drag,  and  “trimmability”  of  the  vehicle  must  be  known  a  priori 
over  the  entire  bight  envelope.  We  present  a  method  that  allows 


for  the  effects  of  a  locked  or  boating  control  effector  to  be 
estimated  over  the  bight  envelope.  This  method  has  the  advantage 
of  including  six  degree-of-freedom  effects  (i.e.,  trimmability)  in 
the  point-mass  models  used  by  trajectory  synthesis  algorithms. 
For  situations  where  there  are  stuck  or  boating  control  effectors, 
the  aerodynamic  data  does  not  change;  therefore,  the  effects 
of  the  failure  can  be  estimated  over  the  entire  bight  envelope. 
We  will  show  how  to  determine  constraints  on  the  trajectory 
using  the  PLCA  approach  that  was  outlined  above.  However,  the 
“preference  vector”  used  in  the  control  allocator  requires  careful 
consideration  along  with  the  specibc  type  of  trajectory  desired 
(e.g.,  computation  of  the  maximum  downrange  vs.  minimum 
downrange  for  the  footprint),  otherwise  adverse  interactions 
between  the  control  allocator  and  the  trajectory  synthesizer  may 
occur. 

I.  Principles  of  Control  Allocation 

Our  interest  in  control  allocation  is  motivated  by  the  de¬ 
velopment  of  re-configurable  control  laws  for  autonomous 
aerospace  vehicles.  One  such  class  of  vehicles  is  the  next 
generation  of  re-usable  launch  vehicles.  The  next  generation 
of  RLVs  will  be  characterized  by  a  minimal  control  effector 
suite  to  facilitate  control  re-configuration  when  the  aircraft 
experiences  either  failures  in  the  flight  control  system  or 
damage  to  the  vehicle.  Typically,  at  the  inner-most  control 
loop,  the  control  system  is  controlling  body-axis  angular  rates 
to  meet  the  outer-loop  guidance  commands.  When  the  number 
of  effectors  exceeds  the  number  of  controlled  variables,  it  is 
quite  common  that  the  desired  commands  can  be  achieved  in 
many  different  ways.  A  control  allocation  algorithm  can  be 
used  to  find  a  set  of  control  effector  positions  that  meet  some 
desired  objective  in  addition  to  delivering  the  desired  angular 
rates.  Additionally,  the  control  effectors  are  subject  to  position 
and/or  rate  limiting  constraints  that  can  be  enforced  by  a  well 
designed  control  allocation  algorithm. 

Control  allocation  can  be  stated  simply  as  follows:  Given 
the  mapping 

ddes  =G(S)  (1) 

subject  to: 

S  G  [fimin  5  fimax\  (2) 

S  G  [fimin  5  fimax\  (3) 
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where  ddes  £  and  5  G  Mn  and  n  >  m,  we  desire  to  find 
the  inverse  mapping  such  that 

5  =  G~\ddes)  (4) 

such  that  the  constraints  6  G  [8min,Smax\  and  S  G 
are  not  violated. 

To  solve  the  control  allocation  problem,  it  is  usually  as¬ 
sumed  that  the  control  moments  produced  by  the  control 
effectors  are  linear  functions  of  the  effector  displacement 

G(S)  =  B6  (5) 

This  has  led  to  much  research  being  devoted  to  various 
methods  for  solving  what  we  will  call  the  linear  control 
allocation  problem.  The  linear  control  allocation  problem  is 
often  re-stated  as  an  equivalent  optimization  problem  (either 
a  linear  or  quadratic  program)  in  order  to  ensure  that  the 
limits  on  the  control  effector  positions  and  rates  are  satisfied 
(see  References  [1],  [2],  [3]).  Non-optimization  approaches  do 
exist,  primarily  those  attributed  to  Durham  [4]. 

The  assumption  that  the  moments  are  linear  functions  of 
the  effectors  is  quite  valid  for  most  control  effectors  about 
their  primary  control  axis,  thus  linear  control  allocation  ap¬ 
proaches  are  successful  under  nominal  flight  conditions.  In 
this  case  small  non-linearities  are  seen  as  disturbances  or 
model  uncertainty  and  the  control  laws  are  robust  enough 
that  system  stability  and  performance  are  not  compromised. 
However,  there  are  situations  where  the  effectors  must  operate 
away  from  the  linear  region,  or  when  it  may  be  beneficial 
to  take  advantage  of  the  non-linear  effects.  In  these  cases, 
we  want  to  include  the  effects  of  the  non-linearities  in  the 
control  allocation.  One  situation  is  when  the  aircraft  has  a 
stuck  effector,  and  the  moment  producing  capability  about  a 
secondary  axis  is  highly  non-linear.  For  example,  for  left-right 
pairs  of  effectors,  there  is  a  quadratic  relationship  between 
the  yawing  moment  and  the  effector  displacement.  In  this 
case,  deflection  of  a  single  effector  on  one  side  of  the  vehicle 
induces  a  yawing  moment  to  that  side  of  the  vehicle  due  to  the 
asymmetric  drag  distribution.  This  yawing  moment  occurs  for 
both  positive  and  negative  deflections  of  the  control  surface. 
In  this  case  a  linear  fit  can  only  be  applied  to  one  side  of  the 
curve,  otherwise  we  have  a  gross  mismodelling  of  the  moment- 
effector  relationship. 

One  approach  that  was  designed  to  account  for  benign 
non-linearities,  and  has  been  successfully  flight-tested,  was 
developed  at  AFRL  by  Oppenheimer  and  Doman  [5].  Their 
approach  added  an  “intercept-correction”  term  that  made  the 
relationship  between  the  controls  and  the  moments  affine: 
B5  +  e  =  ddes •  However,  this  approach  is  still  limiting 
because  there  may  be  slope  reversals  in  the  aerodynamic  data. 
The  possibility  exists  where  the  allocator  could  become  cap¬ 
tured  in  this  non-linear  region,  resulting  in  actuator  saturation. 

To  further  improve  control  allocation  and  to  maximize  the 
performance  of  the  effector  suite  on  the  aircraft  requires  taking 
advantage  of  the  non-linearities  present  in  the  data.  Non-linear 
programming  approaches  were  not  of  interest  due  to  the  fact 
that  polynomials  on  the  order  of  5  or  higher  were  typically 


needed  to  accurately  fit  the  aerodynamic  data,  thus  resulting 
in  a  rather  difficult  set  of  equations  to  solve.  A  more  natural 
approach  is  to  utilize  the  inherent  piecewise  linear  nature  of 
the  aerodynamic  data.  This  is  the  approach  that  we  pursued 
in  Reference  [6],  and  is  outlined  here  for  completeness. 

II.  Piecewise  Linear  Programming 

Piecewise  linear  programming  is  an  optimization  method 
that  allows  one  to  approximate  non-linear  programming  prob¬ 
lems  that  are  comprised  of  separable  functions.  To  solve 
the  resulting  approximation  problem,  a  linear  program  can 
be  formulated  and  then  solved  using  a  modified  simplex 
method  [7].  A  second  option  is  to  formulate  the  non-linear 
program  as  a  mixed-integer  linear  program  [8] . 

In  terms  of  control  allocation,  the  restriction  of  approxi¬ 
mating  separable  functions  by  piecewise  linear  functions  may 
appear  to  be  overly  restrictive.  For  most  aircraft,  the  control  in¬ 
duced  moments  can  be  considered  as  separable  since,  in  many 
cases,  there  are  no  significant  aerodynamic  interactions  among 
the  control  effectors.  In  some  instances  the  cross-coupling  of 
control  effectors  cannot  be  neglected,  such  as  when  control 
effectors  are  located  downstream  of  other  surfaces  or  when 
two  effectors  are  located  adjacent  to  one  another. 

For  purposes  of  illustration,  we  will  approximate  a  single¬ 
valued  function,  f(x),  by  its  piecewise  linear  approximation 
and  show  how  to  formulate  the  minimization  of  /(#),  x  G 
[a,  b\  as  a  piecewise  linear  program.  The  approach  given  below 
for  a  single  variable  function  can  be  generalized  for  multi- 
variable,  separable  functions  rather  easily.  Furthermore,  we 
are  not  restricted  to  only  approximating  the  objective  function 
by  a  piecewise  linear  approximation  since  it  is  also  possible 
to  consider  piecewise  linear  approximations  of  the  constraints, 
if  they  are  separable,  within  the  same  framework.  A  detailed 
discussion  can  be  found  in  Reklaitis,  et.al.  [7]. 

Without  loss  of  generality,  we  begin  by  considering  a 
single  variable  function,  f(x),  defined  on  an  interval,  [a,  b\. 
Begin  by  defining  a  grid  of  K  points  spaced  on  the  interval 
[a,  b]  and  denote  these  points  as  x^k\  k  =  1, . . . ,  K  where 
a  =  x(1>}  <  x <  •  •  •  <  x <  •  •  •  <  x(K^  =  b.  Note 
that  we  are  not  restricted  to  a  uniform  spacing  of  the  x^k\ 
Furthermore,  let  denote  the  value  of  f(x^).  A  piecewise 
linear  approximation  of  f{x)  can  then  be  constructed  by 
connecting  (x^k\  f^)  and  /(fe+1))  with  a  straight 

line  as  shown  in  Figure  1 .  The  equation  of  the  line  connecting 
the  points  (x^k\  fW)  and  (x(/c+1) , is  given  by 

/>)  =f(k)  +  ^fc+i|  _  ~x (II  (®  -  x(k) )  (6) 

where  x  G  [x^ ,  There  will  be  K  —  1  such  equations, 

one  for  each  subinterval.  Observe  that  on  a  given  subinterval, 
x  can  be  written  as 

x  =  AWx(/c)+A(W)x(/c+1)  (7) 

where  >0  and  A^+1)  >  0.  The  A^  are  normalized  such 
that 

A(fe)  +  A(/c+1)  =  1  (8) 
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It  can  then  be  shown  that  Equation  6  can  be  written  as 


/  0)  =  A  (/c)/(/c)  +  A(/c+1)/(/c+1) 

Therefore,  in  the  interval  [x^\  x^],  a  given  value  ol 
the  approximate  value  f{x )  can  be  determined  by  asf 
appropriate  values  to  A^  and  A^+1^  that  correspond 
subinterval  in  which  x  lies.  Since  x  can  only  be  define 
single  subinterval,  all  the  A^  which  are  not  associate 
that  particular  interval  all  must  be  equal  to  zero.  As  a 
we  can  express  Equations  7  and  9  as 

x  =  ^\(k)x(k) 
k=  1 

f(x)  =  j2*(k)f(k) 

k=l 


subject  to  the  following  conditions: 


K 


=  1 

(12) 

k= 1 

\<k)  >0,  k  =  l,...,K 

(13) 

a6)aO)=0  if  j  >*  +  !;*  =  !,. 

(14) 

Equation  14  is  necessary  to  ensure  that  only  points  lying  on 
piecewise  linear  segments  that  connect  adjacent  breakpoints 
are  considered  as  part  of  the  approximating  function.  For 
example,  given  a  value  of  x,  no  more  than  two  of  the  A  ^ ’s 
are  allowed  to  be  non-zero  and  positive,  and  the  two  A  ^ ’s 
also  must  be  adjacent.  If  we  consider  a  value  of  x  where 
A^3)  and  A^4)  are  positive,  with  A^1)  =  A^2^  =  0  and 
\(k)  =  0,  k  =  5, . . . ,  K,  then  the  value  of  f(x)  lies  on  the 
approximating  function  between  x ^  and  x^ .  On  the  other 
hand,  if  A^4^  >0  was  to  be  replaced  by  A^6^  >  0,  and  all  other 
=  0,  then  the  line  connecting  x ^  and  x^  would  not  be 
part  of  the  approximating  function.  Furthermore,  if  we  chose 
a  value  of  x  such  that  x  =  x^  and  f{pc)  =  f(x),  then  from 
Equation  12,  A^  =  1  and  all  other  values  of  A  =  0.  Lastly,  it 
is  important  to  note  that  one  can  always  obtain  a  more  accurate 
approximation  of  f  {x)  by  increasing  the  number  of  gridpoints; 
however,  this  obviously  increases  the  size  of  problem. 

Given  that  we  now  have  a  piecewise  linear  approximation 
to  f{x)  and  the  additional  constraints  that  result  from  the 
transformation,  we  are  able  to  state  the  Piecewise  Linear 
Program  that  corresponds  to  the  minimization  of  f{x)  on  the 
interval  a  <  x  <  b. 

K 

min  f(x)  =  E  A  (fc)/(fe)  (15) 

k= 1 


Fig.  1.  Piecewise  Linear  Function  Approximation 

subject  to 

y^A(fe)  =  i  (i6) 

k=  1 

A(fe)  >  0  (17) 

Once  the  solution  to  the  piecewise  linear  program  is  obtained, 
one  uses  Equation  10  to  find  the  corresponding  value  of  x  that 
gives  an  approximate  minimum  to  f{x).  Finding  a  solution  to 
a  piecewise  linear  program  requires  an  approach  that  ensures 
that  Equation  14  is  satisfied.  Recall  that  Equation  14  requires 
that  no  more  than  two  adjacent  A^’s  are  allowed  to  be 
non-zero.  Therefore,  to  find  an  optimal  feasible  solution  to 
the  piecewise  linear  program,  one  of  two  approaches  must 
be  taken.  One  approach  is  to  solve  the  problem  using  the 
simplex  method  with  a  restricted  basis  entry  rule.  [7]  A 
second  approach  is  to  formulate  Equation  14  using  binary 
decision  variables  [8]  that  will  constrain  x  to  be  on  only  one 
subinterval.  The  result  will  be  yet  another  increase  in  the  size 
of  the  problem  beyond  what  was  necessary  for  the  piecewise 
linear  approximation.  The  addition  of  the  binary  variables 
transforms  the  piecewise  linear  programming  problem  into  a 
mixed-integer  linear  program  (MILP).  We  will  take  the  latter 
approach  because  it  is  sufficient  for  demonstrating  the  validity 
of  the  approach,  and  also  because  of  the  availability  of  an 
open-source  code  (GNU  Linear  Programming  Kit)  that  solves 
linear  programs  and  mixed-integer  linear  programs. 

A.  Transformation  of  the  Piecewise  Linear  Program  to  a 
Mixed-Integer  Linear  Program 

Begin  by  considering  the  piecewise  linear  approximation 
shown  in  Figure  1.  Note  that  if  there  are  K  breakpoints,  then 
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there  are  K  —  1  linear  segments.  We  assign  a  variable  that 
corresponds  to  the  kth  linear  segment  of  the  piecewise  linear 
approximation  such  that 

y{k)  =  f  1  ^  A(fe)  +  0  a01!  A(fe+l)  +  0,  (18) 

[  0  otherwise 

for  k  =  1, . . . ,  K  —  1.  Next,  we  make  the  observation  that  if 
7^  0  and  ^  0,  then 


A(1)  <  y(1)  (19) 

A(2)  <  yw  (20) 

where  ytl>  =  1.  However,  if  we  are  on  the  segment  where 
A^2)  ^  0  and  A®  ^  0,  such  that  y ^  =  1,  then 

A(2)  <  y(2)  (21) 

A  (3)<y(2).  (22) 

If  we  proceed  in  this  manner,  we  observe  that  the  following 

restrictions  can  be  placed  on  the  A^ 

A(1)  <  y(1),  (23) 

A(fc)  <J/(fc-i)+y(*)!  k  =  2, . . . ,  K  —  1  (24) 

X(K)<y(K- l)_  (25) 


The  rationale  behind  Equation  24  is  as  follows:  the  A <k>  that 
correspond  to  points  that  are  interior  to  the  interval  (i.e.,  they 
are  not  the  endpoints  of  the  interval  on  which  x  is  defined)  can 
be  associated  with  one  of  two  line  segments.  A  particular  A^ 
is  an  endpoint  for  the  line  segment  immediately  preceding  it 
in  addition  to  the  line  segment  that  comes  immediately  after  it. 
Only  one  of  these  two  line  segments  may  be  “active”  at  any 
time;  therefore,  the  right-hand  side  of  Equation  24  is  never 
greater  than  one.  In  addition  to  Equations  23-25,  we  have  an 
additional  constraint  to  ensure  that  only  one  of  the  K  —  1  line 
segments  is  active,  hence  only  one  of  the  y  ^  can  be  equal  to 
one: 

K- 1 

£  y{k)  = 1  (26) 

k=  1 

By  including  Equations  23-26  into  the  piecewise  linear 
program  for  the  arbitrary  single  valued  function  being  min¬ 
imized  in  the  example  above,  we  transform  it  into  a  mixed 
integer  linear  program.  The  transformed  optimization  problem 
is  stated  as  follows: 

K 

mmf(x)  =  J2x(k)fik)  (27) 

k=  1 


subject  to 


K 


£a<*>  =1 

(28) 

k—  1 

\{k)  >  0 

(29) 

VI 

(30) 

X(k)  <y{k~l)  +y(k),  k  =  2}. 

(31) 

1 

VI 

g 

(32) 

K- 1 

£^  =  1 

(33) 

k=  1 

y{k)  e  {0, 1} 

(34) 

By  including  the  additional  constraints  which  are  necessary  to 
complete  the  transformation  of  the  piecewise  linear  program, 
we  have  added  an  additional  K  —  1  decision  variables  to 
the  problem.  This  does  not  include  any  slack  or  surplus 
variables  that  may  be  required  by  the  solver  to  convert 
inequality  constraints  into  equality  constraints.  The  slack  and 
surplus  variables  will  further  increase  the  number  of  decision 
variables.  The  solution  to  the  mixed  integer  linear  program 
is  obtained  by  using  a  branch-and-bound  algorithm.  Technical 
details  on  the  branch-and-bound  algorithm  can  be  found  in 
Bertsimas  [8]. 

III.  Piecewise  Linear  Mixed  Optimization  Control 
Allocation 

The  control  allocation  problem  solved  in  Reference  [6]  used 
two  distinct  performance  indices  that  were  similar  to  those 
used  in  Buffington’s  [2]  multi-branch  approach.  However, 
instead  of  a  linear  relationship  between  the  control  effector 
displacements  and  the  control  moments,  a  non-linear  relation¬ 
ship  was  used.  Let  the  non-linear  vector- valued  function  G(d) 
denote  the  relationship  between  the  control  effector  positions 
and  their  moments.  The  function  G(d)  maps  Mn  into  Mm 
where  n  >  m  and  n  is  the  number  of  control  effectors  and 
m  is  the  number  of  controlled  variables.  Let  ddes  denote  the 
controlled  variables.  In  this  case  the  controlled  variables  are 
the  rolling,  pitching  and  yawing  moments.  Buffington’s  multi¬ 
branch  approach  requires  that  two  optimization  problems  be 
solved.  The  first  optimization  problem  is  called  the  control 
deficiency  branch.  The  objective  of  the  control  deficiency 
branch  is  to  minimize  \\W  a(G(S)—ddes)\\i  subject  to  position 
and  rate  constraints  on  d.  The  value  of  the  performance  index 
for  this  optimization  problem  indicates  whether  or  not  ddes 
is  feasible.  If  feasibility  of  the  control  allocation  problem 
has  been  ascertained,  a  second  optimization  problem  is  then 
solved  that  minimizes  some  secondary  objective.  The  objective 
function  of  this  second  optimization  is  typically  taken  to  be 
\\WU(S  —  Sp)\\i  subject  to  G(S)  =  ddes  and  position  and  rate 
constraints  on  S.  This  is  commonly  referred  to  as  the  control 
sufficiency  branch. 

The  mixed  optimization  problem  that  was  formulated  by 
Bodson  [3]  combines  the  two  branches  of  the  multi -branch 
control  allocation  problem  into  a  single  optimization  problem. 
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A  new  parameter  e  is  introduced  for  the  purpose  of  prioritizing 
either  control  deficiency  or  control  sufficiency.  The  mixed 
optimization  problem  is  stated  as: 

min  J  =\\Wa(G(S)  -  ddes)\\%  +  e|| WU(S  -  8P) ||i  (35) 

subject  to: 

8 min  —  8  —  &max  (36) 

where  Wa  =  diag(wai,  wa2,  wa3, . . . ,  wam)  is  a  weight¬ 
ing  matrix  used  to  prioritize  a  given  control  axis,  G(8) 
is  a  non-linear,  vector- valued  function  that  maps  Mn  to 
Mm,  8  is  an  n  x  1  vector  of  control  effectors,  Wu  = 
diag(rcni,  wu2,  wus, . . . ,  wun)  is  a  weighting  matrix  on  the 
control  effectors,  and  8P  is  an  n  x  1  vector  of  “preferred” 
control  effector  displacements.  Again  it  is  assumed  that  n  > 
m.  For  the  moment  we  will  make  the  assumption  that  G(S)  = 
B8.  The  mixed  optimization  problem  can  then  be  posed 
as  a  linear  programming  problem.  The  corresponding  linear 
program  can  then  be  solved  by  any  readily-available  linear 
programming  software. 

Bodson  [3]  gives  one  possible  transformation  to  a  linear 
programming  problem  for  the  optimization  problem  defined  in 
Equations  35  and  36;  however,  we  selected  a  transformation 
approach  that  can  be  found  in  Bertsimas  [8].  The  transforma¬ 
tion  relies  on  the  observation  that  \x\  is  the  smallest  number 
xs  that  satisfies  both  x  <  xs  and  —  x  <  xs.  As  a  result,  we 
are  able  to  pose  the  mixed  optimization  problem  as  follows: 

min  J  =  Wa8s  +  eWuus  (37) 

subject  to: 


8S  >0 

(38) 

us  >  0 

(39) 

B8  -f  8S  ^  ddes 

(40) 

-B8  +  8 s  ^  ddes 

(41) 

8  +  zls  >  8p 

(42) 

8  us  A  8 p 

(43) 

8min  —  8  ^  8max 

(44) 

The  vectors  8S  and  us  are  vectors  of  “slack  variables”.  The 
reason  for  selecting  this  particular  transformation  as  opposed 
to  the  one  in  Bodson  [3]  is  that  this  formulation  allows  us  to 
easily  implement  the  piecewise  linear  function  approximation. 

To  convert  the  linear  programming  problem  into  a  piecewise 
linear  programming  problem,  we  simply  replace  Bd  above 
with  a  piecewise  linear  representation  of  each  control  moment 
curve,  (i.e.,  Li(Si ),  Mi(8i),  and  Ni(Si)).  We  choose  a  set  of 


breakpoints  for  each  Si,  i  =  1, . . . ,  n  such  that: 


Ki 


(45) 

k=  1 

K 

(46) 

k=  1 

Ap)Af)=0,  if  k  >  j  +  1,  j  =  1,  ,. 

.,Ki-l  (47) 

(k) 

where  A)  '  is  a  non-negative  interpolating  coefficient  cor¬ 
responding  to  the  ith  control  effector  at  breakpoint  k,  and 
Ki  denotes  the  number  of  breakpoints  for  the  ith  control 
effector.  Equation  47  is  necessary  in  order  ensure  that  Si  is 

(k) 

approximated  by  no  more  than  two  adjacent  values  of  A  •  \  If 

(k) 

Si  falls  at  a  breakpoint,  there  will  only  be  one  value  of  A)  ; 
that  is  non-zero  and  Equation  47  is  still  valid. 

For  inner-loop  flight  control,  we  assume  that  our  objective 
is  to  find  a  set  of  control  deflections  S  such  that  a  specified 
rolling  moment  L ,  pitching  moment,  M,  and  yawing  moment 
N  are  simultaneously  achieved.  The  piecewise  linear  approxi¬ 
mations  for  the  control  moments  as  a  function  of  S  t  are  written 
as 

Li=L[bi)=Yj\f)Lf)  (48) 

k=l 

Ki 

Mi  =  M (Si)  =J2  V'fe) Mik)  (49) 

k= 1 

Ki 

Ni  =  N(Si)  =  ^k)N^k)  (50) 

k= 1 


where  L^k\  M^k\  and  N-k^  are  the  values  of  the  rolling, 
pitching,  and  yawing  moment  curves  evaluated  at  the  kth 
breakpoint  for  the  ith  control  effector.  We  are  now  able  to 
replace  BS  with  BA  where 


rx(i} 

■  ■  L  ik)  • 

f(Kn 

•  •  X/n 

B  = 

m[2)  . 

. .  M^k)  . 

..  MnK‘ 

n[2)  . 

■  ^k)  . 

..  NiK' 

The  vector  A  is  of  length  Y^=i  anc*  B  is  a  matrix  of  size 
nc  x  Y^i= i  Ki  where  nc  is  the  number  of  controlled  variables. 
In  the  piecewise  linear  optimization  problem,  the  constraints 
8 min  <  8  <  8 max  are  replaced  by  A-^  >  0.  The  upper  and 
lower  bounds  on  8  are  accounted  for  in  the  selection  of  the 
breakpoints  for  each  Si.  Once  we  obtain  an  optimal  solution 
to  the  problem,  we  compute  each  Si  using  Equation  45.  It  is 
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also  necessary  to  include  in  the  problem  the  n  constraints  that 

correspond  to  Equation  46. 

The  resulting  optimization  problem  is 

min  J  =Wa8s  A  eWuus 

(53) 

subject  to: 

8S  >0 

(54) 

us  >  0 

(55) 

BA.  A  8S  A  ddes 

(56) 

BA  A  8S  A  ddes 

(57) 

Ki 

jk)8^k)  A  uSii  >  8Pii,  i  = 

(58) 

k=  1 

Ki 

-  + us,%  >  i  =  i,  . . . , n 

(59) 

k=  1 

x\k^  >0,  i  =  1, . . . ,  n,  k 

=  1 

(60) 

Ki 

y2xik) mi>  i  =  v-  - > n 

(61) 

k=  1 

)  =0,  if  k  A  j  +  1,  j  = 

1 . Ki~l 

(62) 

The  constraint  given  by  Equation  62  is  enforced  using  bi¬ 
nary  variables  and  a  set  of  explicit  constraints  on  the  segments 
of  each  piecewise  linear  control  moment  as  described  above. 

The  resulting  optimization  problem  is  then  a 
linear  program.  The  additional  constraints  are 

mixed-integer 

T— 1 

II 

•<s> 

VI 

"Af 

(63) 

,(k)  .  (k- 1)  .  (fc)  .  .  , 

A \  <y\  +  y\  ,  1  =  1 , . . . ,  to,  k  = 

2,  •  •  • ,  Ki  —  1 

(64) 

s»S*-‘\  <  =  1 . ™ 

(65) 

Ki- 1 

=  <-l . m 

(66) 

k= 1 

Vik)  e{o,i} 

(67) 

IV.  Applications 

There  are  two  important  applications  for  control  allocation 
algorithms  that  take  into  account  the  non-linear  relationship 
between  the  control  effectors  and  the  control  moments.  Our 
original  motivation  for  the  use  of  non-linear  control  allocation 
was  adaptive/reconfigurable  control  systems  as  a  means  to 
maximize  the  performance  by  optimally  selecting  the  control 
deflections  to  meet  the  commanded  body-axis  rotational  rates. 
As  will  be  discussed  below,  the  benefits  of  non-linear  control 
allocation  occur  when  there  are  effector  failures  that  require 
the  remaining  control  effectors  to  operate  in  their  non-linear 
regions.  In  this  section,  we  will  compare  the  performance  of 
the  piecewise  linear  control  allocation  approach  to  a  linear 
allocation  method  in  a  simulation  of  a  re-entry  vehicle  that 


Fig.  2.  Block  Diagram  of  Inner-loop  Dynamic  Inversion  Control  Law 


uses  a  dynamic  inversion  control  law  for  inner-loop  control. 
A  second  application  of  non-linear  control  allocation  is  the 
determination  of  constraints  for  use  with  trajectory  reshaping 
algorithms  [9].  Here  it  is  necessary  to  be  able  to  accurately 
determine  the  range  of  angles-of-attack  over  which  the  vehicle 
can  be  trimmed  in  the  presence  of  control  effector  failures.  An 
additional  benefit  of  solving  the  constraint  estimation  problem 
is  that  trim  lift  and  drag  effects  can  be  included  in  the  tra¬ 
jectory  reshaping  problem  with  little  additional  computational 
burden.  The  trim  maps  and  the  lift  and  drag  perturbations  can 
re  determined  off-line  and  stored  on-board  in  the  guidance 
computer  to  facilitate  real-time  trajectory  re- shaping. 

A.  Dynamic  Inversion  Flight  Control 

Dynamic  inversion  controllers  attempt  to  cancel  and  replace 
the  dynamics  of  the  plant  being  controlled  with  a  set  of  desired 
dynamics.  If  the  fidelity  of  the  on-board  reference  model  is 
high  enough,  then  the  dynamic  inversion  control  law  results 
in  a  closed-loop  system  that  behaves  like  a  decoupled  bank  of 
integrators.  In  the  context  of  flight  control,  a  common  objective 
of  a  dynamic  inversion  control  law  is  to  provide  good  body- 
axis  angular  rate  tracking. 

It  is  assumed  that  a  pilot  or  an  outer-loop  guidance  system 
generates  body-axis  angular  velocity  commands  Pc,  Qc ,  Rc. 
The  inner-loop  dynamic  inversion  control  law  is  designed  such 
that  the  aircraft  tracks  these  angular  velocity  commands  (see 
Figure  2).  The  rotational  dynamics  for  an  aircraft  can  be 
written  as: 

I  Co  =  Gb  —  Cl ?  x  loo  (68) 

where  oo  =  [P  Q  R]T ,  I  is  the  moment-of-inertia  tensor, 
and  G  b  are  the  moments  acting  on  the  vehicle.  We  can  express 
Gb  as  a  sum  that  includes  moments  that  are  due  to  the  wing- 
body  aerodynamics  and  propulsion  system,  which  we  will 
collectively  refer  to  as  the  base  moments,  and  moments  due 
to  the  control  effectors: 


Gb  —  GbaseiiO,  P)  -\-G(P ,  (5)  — 

~L~ 

M 

A 

-L(sy 

M(5) 

N 

base 

N(6) 

where  Gbase(oo,  P)  is  the  moment  generated  by  the  base 
engine-aerodynamic  system  and  G(P,  8)  is  the  sum  of  the 
moments  produced  by  the  control  effectors.  The  parameter 
vector  P  denotes  measurable  or  estimable  quantities  that 
influence  the  body  angular  accelerations  and  includes  variables 
such  as  Mach  number,  angle  of  attack,  sideslip  angle  and 
vehicle  mass  properties  such  as  moments  of  inertia.  Thus  we 
define 

f(u>,  P)  =  Gbase  (w,  P)  -  u  x  Iuj  (70) 
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The  model  used  for  the  design  of  the  dynamic  inversion  control 
law  then  becomes: 

Iu  =  f(u,P)±G{P,8)  (71) 

and  our  objective  is  to  find  a  control  law  that  provides  direct 
control  over  cJ  such  that  cJ  =  elides,  i.e., 

Iu>des  =  f(u>,P)  +  G(P,6 )  (72) 

therefore,  the  inverse  control  must  satisfy: 

I*des-f(t*,P)  =  G(P,6)  (73) 

Since  there  are  more  control  effectors  than  controlled  vari¬ 
ables,  a  control  allocation  algorithm  must  be  used  to  obtain 
a  solution.  Solution  of  this  control  allocation  problem  will  be 
discussed  in  detail  in  a  later  section.  Equation  73  states  that  the 
control  effectors  are  to  be  used  to  correct  for  the  difference 
between  the  desired  accelerations  and  the  accelerations  due 
only  to  the  base  moments. 

B.  Simulation  of  Dynamic  Inversion  Control  Law 

The  mixed-integer  linear  programs  for  the  mixed- 
optimization  control  allocation  discussed  above  were  imple¬ 
mented  in  a  Simulink  simulation  of  a  re-entry  vehicle.  This 
particular  vehicle  has  six  control  surfaces:  left  and  right 
rudders,  left  and  right  flaperons,  a  body  flap,  and  a  speed 
brake.  The  simulation  models  the  descent,  final  approach,  and 
touchdown  of  the  vehicle. 

The  performance  of  the  piecewise  linear  approach  is  com¬ 
pared  to  that  of  a  linear  control  allocation  method.  The  linear 
control  allocator  assumes  that  the  moments  are  linear  functions 
of  the  effectors.  The  slope  of  the  control  moment  curve  is 
calculated  with  respect  to  the  current  control  effector  position 
using  a  forward  difference  approximation.  An  intercept  correc¬ 
tion  [5]  term  is  then  applied  to  account  for  mild  non-linearities 
in  the  aerodynamic  data. 

The  results  that  follow  give  the  closed-loop  vehicle  per¬ 
formance  when  there  are  two  failures  injected  into  the  flight 
control  system  at  different  times  during  the  approach  and 
landing  phases.  It  is  assumed  that  there  is  some  type  of 
fault  detection  capability  on-board  the  aircraft  to  identify 
the  failures.  The  failure  information  is  immediately  passed 
to  the  control  allocation  algorithm  in  order  to  facilitate  re¬ 
configuration  of  the  vehicle’s  effectors.  The  aircraft’s  tra¬ 
jectory  begins  at  an  altitude  of  about  15,000  ft  above  the 
runway  and  4  miles  downrange  from  the  runway  threshold. 
The  first  failure  occurs  30  s  into  the  simulation,  and  involves 
the  body  flap  being  locked  at  —5  deg.  This  failure  contributes 
a  constant  pitching  moment  to  the  aircraft.  A  second  failure, 
where  the  right  rudder  becomes  locked  at  1  deg,  occurs  at 
40  s.  This  particular  failure  adds  not  only  a  pitching  moment 
to  the  aircraft,  but  also  rolling  and  yawing  moments.  This 
particular  failure  combination  was  chosen  because  it  requires 
the  flaperons  to  operate  in  a  highly  non-linear  region  of  the 
control  moment  curve.  After  the  failures  are  introduced,  the 
aircraft  tries  to  follow  the  nominal  approach  trajectory  to  the 
runway  threshold.  The  aircraft  extends  the  landing  gear  at 


about  68  s  and  flares  immediately  before  touchdown.  The 
simulation  ends  at  touchdown  when  the  weight-on-wheels 
switch  is  triggered. 

For  the  control  sufficiency  branch  of  the  control  allocator, 
a  “preferred”  control  position,  Sp,  is  required.  There  are 
several  different  objectives  that  may  be  used  to  determine 
Sp.  These  include,  but  are  not  limited  to,  minimum  control 
deflection  (Sp  =  0),  minimum  2-norm  of  deflection,  and  null- 
space  injection  [2],  [10].  The  preference  vector  used  in  this 
simulation  is  the  minimum  2-norm  of  the  control  surface 
deflection.  This  particular  Sp  minimizes  STWS  subject  to 
BS  =  ddes  where  W  is  a  positive  definite  weighting  matrix. 
For  our  results,  we  take  the  weighting  matrix,  W,  as  the 
identity  matrix.  The  corresponding  solution  to  this  problem 
is  then  Sp  =  BT (BBT)~1ddes-  Note  that  this  particular 
preference  vector  has  the  advantage  of  facilitating  robustness 
analysis  with  the  control  allocator  in  the  loop  since  the  control 
allocator  can  be  represented  in  closed-form  for  local  linear 
analysis  .  For  the  piecewise  linear  control  allocation,  we  found 
that  it  was  sufficient  to  compute  the  right  pseudo-inverse 
solution,  Sp,  with  a  B  matrix  that  uses  local  slopes  of  the 
control  moments  at  the  last  control  surface  position. 

We  will  measure  the  performance  of  the  two  control  alloca¬ 
tors  by  their  ability  to  produce  deflections,  that  when  applied 
to  the  non-linear  aerodynamic  database,  produce  the  desired 
moments  about  each  axis.  This  metric  is  an  indication  of  the 
error  that  results  from  the  selection  of  a  particular  model  in 
the  control  allocation  algorithm. 

1 )  Simulation  Results:  The  results  for  the  piecewise  linear 
control  allocator  as  compared  to  a  linear  control  allocator  with 
intercept  correction  are  given  in  Figures  3  and  4.  Figure  3 
shows  the  base  10  logarithm  of  || ddes  ~  G(P,S)  || 2,  where 
G(P,  8)  is  the  moment  that  is  applied  to  the  vehicle  when 
given  the  control  deflections  returned  by  the  control  allocator. 
It  is  evident  that  the  piecewise  linear  control  allocator  returns 
control  surface  deflections  that  produce  the  desired  moments. 
On  the  other  hand,  the  linear  control  allocator  has  a  significant 
error.  Note  that  at  the  40s  mark,  when  the  second  failure 
is  introduced,  both  control  allocators  indicate  that  there  is  a 
moment  deficiency  due  to  control  effector  saturation.  Beyond 
60s  the  performance  of  the  piecewise  linear  control  allocator 
improves  once  the  effectors  are  no  longer  saturated.  The  poor 
performance  of  the  linear  control  allocator  is  primarily  due 
to  the  modelling  errors  inherent  in  the  linear  approximation 
of  the  control  moment  curves.  The  control  surface  commands 
from  each  control  allocator  are  shown  in  Figure  4.  We  see 
that  after  the  second  failure  is  injected  into  the  simulation  that 
the  flaperon  deflections  saturate  at  their  upper  limit.  Note  that 
the  flap  and  speedbrake  commands  for  the  piecewise  linear 
allocator  oscillate  after  about  65  sec.  This  appears  to  be  due 
to  an  oscillatory  ddes  and  Sp  that  result  from  body-axis  rate 
loop  closures.  For  the  linear  control  allocator,  we  see  that 
the  flaps,  left  rudder,  and  speedbrake  exhibit  large  amplitude 
oscillations.  Note  that  for  this  failure  case,  the  vehicle  that 
flew  with  the  piecewise  linear  control  allocator  maintained 
controlled  flight. 
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Fig.  3.  Difference  Between  Desired  and  Applied  Moments 
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ranges  of  angle-of-attack,  for  a  given  Mach  number,  at  which 
the  aircraft  can  be  trimmed.  Also,  it  is  important  to  be  able 
to  estimate  critical  parameters,  such  as  the  maximum  lift- 
to-drag  ratio,  that  will  impact  the  optimal  trajectory  [9].  It 
is  not  desirable  to  use  a  linear  control  allocation  approach 
to  determine  whether  or  not  the  vehicle  can  be  trimmed. 
The  reason  is  that  a  linear  control  allocator  would  require 
a  single  line  of  best  fit  to  approximate  the  moment-deflection 
curve  at  a  fixed  flight  condition,  thereby  introducing  significant 
modelling  errors.  These  modelling  errors  would  then  result 
in  an  inaccurate  determination  of  the  range  of  trimmable 
angles-of-attack.  On  the  other  hand,  using  a  non-linear  control 
allocation  algorithm  allows  one  to  more  accurately  estimate  the 
range  of  angles-of-attack  over  which  one  may  trim  the  aircraft 
as  well  as  the  critical  parameters  that  affect  the  trajectory,  such 
as  trim  lift  and  trim  drag. 

To  solve  for  the  range  of  angles-of-attack  over  which  the 
vehicle  can  be  trimmed,  we  want  to  determine  where  in  the 
flight  envelope  there  exists  enough  available  control  power  to 
null  the  wing-body  rolling,  pitching,  and  yawing  moments. 
The  procedure  for  this  is  outlined  below  in  Algorithm  1.  If 
we  consider  Equation  73,  and  set  u>  =  Co  =  0,  then  we  have 
—  /(0,P)  =  G(P,5).  The  desired  moment  vector  will  then 
be  defined  as  ddes  =  —  /(0,  P).  Thus,  we  have 


(74) 


Fig.  4.  Control  Effector  Time  History 
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It  should  be  noted  that  both  the  linear  and  non-linear  control 
allocation  approaches  were  solved  using  the  GNU  Linear 
Programming  Kit’s  linear  programming  simplex  library  for 
the  former  and  the  branch-and-bound  solver  for  the  latter. 
The  simulation  with  the  non-linear  control  allocator  and  the 
branch-and-bound  solver  ran  extremely  slow  when  compared 
to  the  simulation  with  the  linear  program  solver  used  for 
the  linear  control  allocator.  It  was  observed  that  it  took 
approximately  8.5  times  longer  on  average  to  find  a  solution 
to  the  MILP  control  allocation  problem  as  compared  to  the 
linear  programming  approach.  This  computational  burden  may 
limit  one’s  ability  to  utilize  the  MILP  approach  in  a  real-time, 
digital  flight  control  system  in  the  near  future.  Performance 
gains  may  be  achieved  by  solving  the  piecewise  linear  control 
allocation  problem  via  the  simplex  method  with  the  restricted 
basis  entry  rules  [7]. 

C.  Application  to  Trajectory  Reshaping 

The  intent  here  is  to  show  a  second  important  application 
of  non-linear  control  allocation.  One  of  the  areas  of  active 
research  regarding  reusable  launch  vehicles  is  the  on-line 
determination  of  new  feasible  trajectories  following  a  control 
effector  failure.  It  is  desired  to  safely  guide  the  aircraft  from 
the  time  that  the  failure  occurs  to  a  safe  abort  and  to  recover 
the  vehicle  if  possible.  An  important  part  of  the  determination 
of  feasible  trajectories  is  having  accurate  estimates  of  the 


where  Lbase  is  the  wing-body  rolling  moment,  Mbase  is  the 
wing-body  pitching  moment,  and  Nbase  is  the  wing-body 
yawing  moment.  Typically  symmetric  flight  (i.e.,  zero  sideslip 
angle)  is  assumed;  therefore,  the  side  force  and  Lbase  and 
Nbase  are  taken  to  be  zero.  Note  the  presence  of  the  minus 
sign  above  in  Equation  74.  The  minus  sign  is  present  because 
we  want  to  counteract  the  wing-body  moment  with  the  control 
effectors  to  achieve  rotational  equilibrium;  therefore,  we  desire 
the  moment  produced  by  the  effectors  to  be  of  opposite  sign 
to  that  of  the  base  moments.  More  details  on  the  algorithm 
can  be  found  in  Reference  [11]. 

Shown  in  Ligure  5  is  an  example  of  the  moment  deficiency 
for  the  X-33  vehicle  model  used  by  Shaffer  [12].  The  X-33 
has  a  suite  of  8  effectors:  right  and  left  body  flaps,  right  and 
left  rudders,  and  right  and  left  inboard  and  outboard  elevons. 
In  Ligure  5  the  contours  of  the  trim  deficiency  as  a  function 
of  Mach  and  a  are  given.  The  failure  that  is  considered  is  a 
double  body-flap  stuck  at  26  deg  trailing  edge  down.  A  non¬ 
zero  value  of  the  control  deficiency,  \\ddes  —  G(S)  ||2,  indicates 
that  vehicle  cannot  be  trimmed  at  that  particular  Mach  number 
and  angle-of-attack.  Notice  that  the  region  over  which  the 
aircraft  can  be  trimmed  is  significantly  reduced.  Only  below 
bout  Mach  0.25  can  the  vehicle  be  trimmed  at  any  angle-of- 
attack.  Ligures  6  and  7  demonstrate  the  change  in  lift  and  drag 
over  the  flight  envelope.  Note  that  a  negative  AC d  indicates 
an  increase  in  drag  and  a  negative  AC  l  indicates  an  increase 
in  lift  compared  to  the  un-failed  case. 

Of  particular  importance  when  computing  the  trim  map 
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Algorithm  1  Trim  Map  Generation 

Define  a  grid  for  Mach  number  and  angle-of-attack  and 
discretize  such  that  M  =  {Mi, . . . ,  Mi , . . . ,  Mm}  and 

OL  =  .  OLj ,  .  .  .  ,  C^tv} 

for  i  =  1  to  M  do 
for  j  =  1  to  AT  do 

Determine  the  wing-body  moments,  ddes  = 
0 

M-base  (Mj,  OLj  ) 

L  o 

Pick  the  preference  vector, 

Solve  the  control  allocation  problem  for  £  given 
and 

Compute  the  control  deficiency 
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Fig.  6.  Lift  Coefficient  Perturbation  Due  to  Body  Flaps  Failed  at  26  deg 
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Change  in  Drag  Coefficient 


Compute  the  lift  and  drag  at  ( Mi ,  ay)  using  5  obtained 
above. 

end  for 
end  for 


Pitch  Deficiency  Contour 


-10  0  10  20  30  40  50 

Alpha  (deg) 


Fig.  5.  Moment  Deficiency  for  a  26  deg  Body  Flap  Failure 


is  the  selection  of  the  preference  vector,  Sp.  The  moment 
deficiency  is  independent  of  the  preference  vector  since  it 
is  determined  solely  by  the  control  deficiency  branch  of  the 
cost  function:  \\ddes  ~  G(5,P)\\i.  At  the  trimmable  Mach- 
a  points,  there  may  exist  non-unique  solutions  to  the  control 
allocation  problem,  thus  the  trim  effects  on  the  lift  and  drag 
are  effected  by  the  deflection  of  the  remaining  effectors. 
Computing  the  lift  and  drag  effects  allows  us  to  capture  6 
degree-of-freedom  effects  in  the  reduced  order  models  used 
for  trajectory  synthesis.  Figures  8  and  9  show  the  effects 
of  the  preference  vector  on  the  lift  and  drag  coefficients  for 
the  nominal,  unfailed  case.  The  two  preference  vectors  being 
compared  are  a  minimum  drag  preference  vector,  upref  =  0, 
and  a  maximum  drag  case  where  upref  was  set  at  the  upper 
limit  for  all  effectors.  Figures  8  and  9  indicate  the  level  of 
sensitivity  of  the  lift  and  drag  coefficient  to  the  preference 
vector. 

The  objective  that  one  is  trying  to  obtain  when  performing 
trajectory  synthesis  must  be  kept  in  mind  when  selecting  the 


Fig.  7.  Drag  Coefficient  Perturbation  Due  to  Body  Flaps  Failed  at  26  deg 
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Fig.  8.  Change  in  Lift  Coefficient  Due  to  Preference  Vector  for  No  Failures 


Change  in  Drag  Coefficient 


Fig.  9.  Change  in  Drag  Coefficient  Due  to  Preference  Vector  for  No  Failures 
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preference  vector.  An  obvious  situation  where  an  inappropriate 
choice  of  the  control  preference  vector  in  the  control  allocation 
can  have  a  negative  implication  is  when  determining  the  land¬ 
ing  footprint  [12].  For  example,  selecting  a  preference  vector 
to  minimize  the  drag  on  the  vehicle  when  trying  to  compute 
the  minimum  downrange  will  produce  an  erroneous  solution. 
In  order  to  minimize  downrange,  it  is  desired  to  minimize  total 
L/D  for  the  vehicle.  Therefore,  the  control  allocator  should 
be  set  to  meet  the  moment  demand  while  increasing  the  trim 
drag  and  decreasing  the  trim  lift.  Conversely,  to  maximize 
the  downrange,  one  wants  to  maximize  L/D.  Therefore,  to 
minimize  the  trim  drag,  the  preference  vector  should  be  set 
to  dp  =  0.  Poor  selection  of  the  preference  vector  will  result 
in  the  control  allocation  having  unintended,  adverse  effects  on 
the  trajectory. 

V.  Conclusions 

A  method  was  presented  for  the  solution  of  a  class  of  non¬ 
linear  control  allocation  problems.  Control  allocation  has  his¬ 
torically  been  performed  by  assuming  that  a  linear  relationship 
exists  between  the  control  induced  moments  and  the  control 
effector  displacements.  Since  aerodynamic  data  almost  always 
exhibits  non-linear  behavior,  such  assumptions  can  lead  to 
degraded  performance  or  vehicle  loss  when  secondary  non¬ 
linear  effects  must  be  used  to  control  a  vehicle,  particularly 
after  control  effector  failures  have  occurred.  Aerodynamic 
databases  are  usually  discrete-valued  and  are  almost  always 
stored  in  multi-dimensional  look-up  tables  where  it  is  assumed 
that  the  data  is  connected  by  piecewise  linear  functions. 

The  approach  that  was  presented  assumes  that  the  control 
effector  moment  data  is  piecewise  linear.  This  assumption 
allows  us  to  cast  the  control  allocation  problem  as  a  piecewise 
linear  program.  In  order  to  solve  the  piecewise  linear  program, 
it  was  re-formulated  as  a  mixed-integer  linear  program  and 
solved  using  a  readily  available  branch-and-bound  algorithm. 
Simulation  showed  that  the  piecewise  linear  programming 
formulation  results  in  improved  tracking  performance  of  the 
desired  moments  when  compared  to  a  more  traditional  control 
allocation  approach  that  uses  the  linear  assumption,  especially 
when  the  aircraft  is  forced  to  operate  with  its  control  ef¬ 
fectors  in  non-linear  portions  of  the  control  moment  curves 
as  a  result  of  control  effector  failures.  The  piecewise  linear 
control  allocator  was  able  to  maintain  control  of  the  aircraft 
and  land  after  the  failures  were  introduced  while  a  control 
allocation  algorithm  that  utilized  a  simple  linear  relationship 
along  with  an  intercept  correction  term  did  not.  The  piecewise 
linear  control  allocator  was  applied  to  the  determination  of 
regions  of  trimmable  angle-of-attack  and  Mach  number  for  the 
purposes  of  trajectory  reshaping  under  failure  conditions.  We 
also  discussed  the  effects  of  the  choice  of  preference  vector  to 
estimate  trim  lift  and  trim  drag  in  conjunction  with  the  reduced 
order  trajectory  optimization  algorithms. 
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